On the Formation of Copper Linear Atomic Suspended Chains 
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We report high resolution transmission electron microscopy and classical molecular dynamics sim- 
ulation results of mechanically stretching copper nanowires conducting to linear atomic suspended 
chains (LACs) formation. In contrast with some previous experimental and theoretical work in 
literature that stated that the formation of LACs for copper should not exist our results showed the 
existence of LAC for the [111], [HO], and [100] crystallographic directions, being thus the sequence 
of most probable occurence. 

PACS numbers: 66.30.Pa,68.37.Lp,68.65.-k,81.07.Lk 



In the last years an enormous amount of theoreti- 
cal and experimental effort have been devoted to the 
study of nanostructures. The advances of experimental 
high spatial resolution techniques or microscopies (high 
resolution transmission electron microscopy (HRTEM), 
scanning tunneling microscopy (STM), and atomic force 
microscopy (AFM)) have revealed new phenomena at 
nanoscale. 

Among nanostructures, metallic nanowires (NWs) are 
of great interest due to the observation of very inter- 
esting physical phenomena (spin filters, quantized con- 
ductance, etc.) and their possible technological applica- 
tion in nanodevices (molecular electronics, nanocontacts, 

etc.) Hiisiiiiiiaiiiiiiiiiiiiiii. 

These structures have displayed quantized conductance 
of n le^jh (Go), where n is an integer number, e is 
the electron charge and, h is Planck's constant 0. The 
ultimate NWs are one-atom thick, when suspended lin- 
ear atomic chains (LACs) are formed by the mechanical 
stretching. 

Most studies of LACs structures have been focused on 
Au, but in principle LACs may also possible to form 
for many other fee metals (Pt,Ag,Pd,etc.), although we 
must emphasize that this issue is still rather controver- 
sial m E 113. In this work, we present a de- 
tailed theoretical (molecular dynamics simulations) and 
experimental (HRTEM) study of mechanical stretching 
of copper NWs, in particular, addressing the formation 
of suspended atom chains. 

Metalhc NWs were produced in situ in the HRTEM 
(JEM-3010 URP 300kV, 0.17 nm point resolution) us- 
ing the methodology proposed by Kondo and Takayanagi 
|21| . Initially holes are opened at several points in a self- 
supported metal film by focusing the microscope elec- 
tron beam (~120 A/cm^). When two holes become very 
close, nanometric constrictions (bridges) are formed be- 
tween them. Then, the microscope beam current derisit 



is reduced to '--^10-30 A/cm for image acquisitio: 
and the nanowires evolve spontaneously, elongate and fi- 
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nally break, sometimes with the appearance of atomic 
suspended chains (LACs). These processes are regis- 
tered using a high sensitive TV camera (Gatan 622SC, 
30 franies/s) and a standard video recorder. 

The polycrystalline Cu thin films (10-30 nm in thick- 
ness) were prepared by thermal evaporation of metal in 
a standard vacuum evaporator (lO"'' mbar). A quartz 
crystal monitor was used to set the evaporation rate of 
metal source and, subsequently, to measure the equiva- 
lent thickness of the film. To prevent oxidation the films 
were sandwiched between two (~ 3 nm thick) amorphous 
carbon layers. Inside the microscope, the carbon layers 
are removed by electron irradiation |l8| |. 

In Fig. 1, we present a typical HRTEM snapshot show- 
ing a LAC composed by two atoms. These represent a 
direct atomic resolution observation of LAC for copper 
(see movie 01) Due to copper relative low atomic 

number and smaller lattice parameter, HRTEM image 
contrast is rather weak, then, the observation of high 
contrast atomic resolution image is much more difficult 
to obtain than for other metals such as Au and Pt [l8j |. 
Experimentally, LAC formation has been observed for 
NWs stretched along the [100], [110], and [111] crystal- 
lographic directions (2J] . 

In order to address the dynamical aspects of the copper 
LAC formation we have also carried out molecular dy- 
namics simulation using a methodology we have recently 
developed and, that has been proved to be very effective 
to address different aspects of metal NW elongation and 
rupture |U ■ Our methodology is based on tight- 

binding molecular dynamics (TB-MD) techniques j2g 
using second-moment approximation (SMA) |27| with 
a small set of adjustable parameters. Details of the 
methodology were published elsewhere |.25i] . In this work, 
the physical stretching that generates the NWs is sim- 
ulated through structural modification of the wire size 
(increasing the distance between the outmost layers, for 
instance). We have considered different elongation rate 
variations, different temperatures, and initial velocity 
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FIG. 1: Snapshot from HRTEM showing LAC with two hang- 
ing and between two grains with different crystallographic ori- 
entation along the elongation direction, [110] on left and [100] 
on right side. Note that the atomic planes of the right apex 
show the expected d2oo ~ 0.18 nm lattice fringes spacing. 



distributions. 

The initial NW configurations consisted of NWs ori- 
ented along [111] (hereafter noted [111] NWs), [110], and 
[100] crystallographic directions with the number of par- 
ticles varying between 450-900 copper atoms (typically 
5-9 atoms in square basis length and with 10-12 layers 
along the elongation direction). In order to mimic the 
bulk structure influence upon NW behavior, we have used 
buffer atom layers as contacts depending on the crystal- 
lographic orientations, two for [110] and [100], and three 
for [111] in order to preserve the AB(ABC) stacking layer 
patterns. These buffer layers are frozen and constrained 
to move only along the stretching directions, while the 
remaining atoms are free to move into any directions. 
We have used temperature range of 300-350K, pulling 
velocities between 0.1-1.5 m/s, time step of 1.0 fs and, 
atomic interactions up to fifth neighbors. A total of 600 
MD simulations were performed (200 for each one of the 
three analyzed crystallographic directions). In Fig. 2 we 
present some examples of snapshots from MD simulations 
for different crystallographic directions. 

Our results have showed the possibility of LAC for- 
mation in the final stages of copper NW stretching for 
all crystallographic studied directions. The number of 
atoms in the LAC are dependent, as expected, on the 
pulling velocity, temperature and initial velocity distri- 
bution. As a general trend we have observed that at 
lower pulling velocity (~ O.lm/s) the system has enough 
time to relax and distribute the internal stress. Then, 
a better structural order and faceting of the apexes is 
obtained and, LACs can be well formed as a linear ar- 
rangement of hanging atom between to atomically sharp 
tips. At higher puUing velocity, (~ 1.5m/s), it is possible 
to generate longer LACs, but with a poor apex faceting 
and, an merely approximative linear distribution of sus- 
pended atoms along the chain. The forces acting on the 
LACs atoms just before rupture are about ~1.5nN, what 
is in excellent agreement with experimentally observed 
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FIG. 2: Snapshots from MD simulations of Cu NW elongation 
along (a) [111], (b) [110] and, (c) [100] crystallographic direc- 
tions, pulling velocity at l.Om/s, temperature 300K. LAC for- 
mation can be clearly identified for the [111] NW in (a). Note 
that the wires in (b) and (c) exhibit pentagonal structures 
(see text for discussions) 



TABLE I: Statistics of LACs formation along crystallographic 
orientations. 



Crystal Direction 


[100] 


[110] 


[111] 


Statistics (%) 


1 (2/200) 


55 (110/200) 


65 (130/200) 



forces for gold NWs 12^ ]. 

Table 1 summarizes the statistics of LAC formation de- 
rived form the MD simulations. Although the formation 
of suspended chains has been observed for all directions 
they occur predominantly for [111] and [110] NWs (Figs. 
2a, b). Although Cu, Au and Ag share many common 
structural features, the NWs exhibit clearly differenti- 
ated patterns regarding LAC occurence P, 0] . Experi- 
mentally an theoretically, it has been observed that Au 
chains are generated with decreasing occurrence rate for 
[111], [100] and [110] axes [llli], while Ag show a com- 
pletely different pattern ([100], [110] and, [111], in de- 
creasing order 19]). Our simulations indicate that, Cu 
LACs show behavior which is distict from Au and Ag 
([111], [110] and, [100] in decreasing order). As men- 
tioned above due to intrinsic difficulties associated with 
copper HRTEM experiments, altough we do not have 
enough experimental data to build a sound statistics, for 
the available experimental data the general trends are in 
agreement with our molecular dynamics results. 

Non-crystallographic icosahedric structures have been 
observed for wires along the [110] and [100] directions 
(Fig. 2b and 2c respectively). These kind of atomic 
arrangement has been recently associated with copper 
[110] NW in order to explain electric transport experi- 
ments displaying a 4.5Go quantum conductance plateaus 
[i^j ] . As a co nseq uence of the modified geometrical Wulff 
construction [23], it can expected that [110] NWs form 
icosahedric structures with five fold axes symmetry |23l ]. 
IN order to explain why pentagonal wires may be formed 
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for Cu [100] NWs, we must also consider that the struc- 
tural rearrangement in the apexes, may sometimes gen- 
erate quite asymmetric junctions. This may lead to the 
formation of NWs, which may be tilted in relation to the 
original direction; in fact, the NW in Fig. 2c is tilted by 
~35 degrees from the [100] axis. Rodrigues et al have 
already revealed that the atomic arrangement of atomic- 
size NWs adjusts such that one of the main gold zone 
axes lies approximately parallel to the elongation direc- 
tion, independent of the apex crystal orientation, and 
this structural adjustment occurs mainly concentrated 
around the narrowest neck region. Then, it is not sur- 
prising that a Cu [100] NW tilted by ^^35 degress (shown 
in Fig. 2C) could adopt a [110]-like atomic structure 
(icosahedric), what in fact means a minor 10 degrees dis- 
orientation from a [110] axis. This is quite well repro- 
duced in Fig. 2c, and it is a good indication that the 
simulations are confirming expected structures. 

As discussed above, dynamical features of the neck 
or apex atomic reorganization influence LAC formation 
and, then structural/morphological aspects are very im- 
portant. In this sense the size of crystalline grains (or 
initial configuration for a simulation) might play an sig- 
nificant role. In order to evaluate the importance of the 
aspect ratio of the initial NWconfiguration, we have car- 
ried out further simulations, where we systematically in- 
creased the lateral dimensions. We have observed that 
after certain limit, we start to observe the appearance of 
parallel columnar NWs; in Fig. 3 we present an example 
of this phenomena. The general trends of each constric- 
tion are basically the same as in the case of only one 
collum. 

The existence of this kind of structures was first pro- 
posed by Marquez and Garcia 30]. Although in our 
HRTEM experiments they are unlike to occur due to 
fabrication procedures (thinning a bridge between two 
growthing holes). They are a real possibility in mechan- 
ically controllable break junction (MCBJ) experiments 

Q. 

In summary our HRTEM and molecular dynamics sim- 
ulatiosn have shown that, in contrast with some previ- 
ous theoretical and experimental predictions, copper sus- 
pended linear atomic chains are possible for [111], [110], 
and [100], being this the sequence of likely ocurrencc. 
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FIG. 3: Sequential snapshots (a), (b) and (c) of MD simula- 
tions showing the occurence of muIti-coUumnar structures in a 
system with 1800 copper atoms with pulling velocity (4.0m/s 
alloy y axis) direction along [111] crystallographic direction, 
temperature 350K. Left side pictures are view along x axis 
and right side pictures along z axis. These multi-columnar 
structures could form LAC one- or two- atom-long. 
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